Optimizing selection based on BLUPs or BLUEs in multiple sets of genotypes differing in their population parameters.

KEY MESSAGE: Selection response in truncation selection across multiple sets of candidates hinges on their post-selection proportions, which can deviate grossly from their initial proportions. For BLUPs, using a uniform threshold for all candidates maximizes the selection response, irrespective of differences in population parameters. Plant breeding programs typically involve multiple families from either the same or different populations, varying in means, genetic variances and prediction accuracy of BLUPs or BLUEs for true genetic values (TGVs) of candidates. We extend the classical breeder's equation for truncation selection from single to multiple sets of genotypes, indicating that the expected overall selection response ( Δ G Tot ) for TGVs depends on the selection response within individual sets and their post-selection proportions. For BLUEs, we show that maximizing Δ G Tot requires thresholds optimally tailored for each set, contingent on their population parameters. For BLUPs, we prove that Δ G Tot is maximized by applying a uniform threshold across all candidates from all sets. We provide explicit formulas for the origin of the selected candidates from different sets and show that their proportions before and after selection can differ substantially, especially for sets with inferior properties and low proportion. We discuss implications of these results for (a) optimum allocation of resources to training and prediction sets and (b) the need to counteract narrowing the genetic variation under genomic selection. For genomic selection of hybrids based on BLUPs of GCA of their parent lines, selecting distinct proportions in the two parent populations can be advantageous, if these differ substantially in the variance and/or prediction accuracy of GCA. Our study sheds light on the complex interplay of selection thresholds and population parameters for the selection response in plant breeding programs, offering insights into the effective resource management and prudent application of genomic selection for improved crop development.

Whole-genome resequencing of Chinese indigenous sheep provides insight into the genetic basis underlying climate adaptation.

BACKGROUND: Chinese indigenous sheep are valuable resources with unique features and characteristics. They are distributed across regions with different climates in mainland China; however, few reports have analyzed the environmental adaptability of sheep based on their genome. We examined the variants and signatures of selection involved in adaptation to extreme humidity, altitude, and temperature conditions in 173 sheep genomes from 41 phenotypically and geographically representative Chinese indigenous sheep breeds to characterize the genetic basis underlying environmental adaptation in these populations. RESULTS: Based on the analysis of population structure, we inferred that Chinese indigenous sheep are divided into four groups: Kazakh (KAZ), Mongolian (MON), Tibetan (TIB), and Yunnan (YUN). We also detected a set of candidate genes that are relevant to adaptation to extreme environmental conditions, such as drought-prone regions (TBXT, TG, and HOXA1), high-altitude regions (DYSF, EPAS1, JAZF1, PDGFD, and NF1) and warm-temperature regions (TSHR, ABCD4, and TEX11). Among all these candidate genes, eight ABCD4, CNTN4, DOCK10, LOC105608545, LOC121816479, SEM3A, SVIL, and TSHR overlap between extreme environmental conditions. The TSHR gene shows a strong signature for positive selection in the warm-temperature group and harbors a single nucleotide polymorphism (SNP) missense mutation located between positions 90,600,001 and 90,650,001 on chromosome 7, which leads to a change in the protein structure of TSHR and influences its stability. CONCLUSIONS: Analysis of the signatures of selection uncovered genes that are likely related to environmental adaptation and a SNP missense mutation in the TSHR gene that affects the protein structure and stability. It also provides information on the evolution of the phylogeographic structure of Chinese indigenous sheep populations. These results provide important genetic resources for future breeding studies and new perspectives on how animals can adapt to climate change.

Genomic evidence for human-mediated introgressive hybridization and selection in the developed breed.

BACKGROUND: The pig (Sus Scrofa) is one of the oldest domesticated livestock species that has undergone extensive improvement through modern breeding. European breeds have advantages in lean meat development and highly-productive body type, whereas Asian breeds possess extraordinary fat deposition and reproductive performance. Consequently, Eurasian breeds have been extensively used to develop modern commercial breeds for fast-growing and high prolificacy. However, limited by the sequencing technology, the genome architecture of some nascent developed breeds and the human-mediated impact on their genomes are still unknown. RESULTS: Through whole-genome analysis of 178 individuals from an Asian locally developed pig breed, Beijing Black pig, and its two ancestors from two different continents, we found the pervasive inconsistent gene trees and species trees across the genome of Beijing Black pig, which suggests its introgressive hybrid origin. Interestingly, we discovered that this developed breed has more genetic relationships with European pigs and an unexpected introgression from Asian pigs to this breed, which indicated that human-mediated introgression could form the porcine genome architecture in a completely different type compared to native introgression. We identified 554 genomic regions occupied 63.30 Mb with signals of introgression from the Asian ancestry to Beijing Black pig, and the genes in these regions enriched in pathways associated with meat quality, fertility, and disease-resistant. Additionally, a proportion of 7.77% of genomic regions were recognized as regions that have been under selection. Moreover, combined with the results of a genome-wide association study for meat quality traits in the 1537 Beijing Black pig population, two important candidate genes related to meat quality traits were identified. DNAJC6 is related to intramuscular fat content and fat deposition, and RUFY4 is related to meat pH and tenderness. CONCLUSIONS: Our research provides insight for analyzing the origins of nascent developed breeds and genome-wide selection remaining in the developed breeds mediated by humans during modern breeding.

Adaptive introgression reveals the genetic basis of a sexually selected syndrome in wall lizards.

The joint expression of particular colors, morphologies, and behaviors is a common feature of adaptation, but the genetic basis for such "phenotypic syndromes" remains poorly understood. Here, we identified a complex genetic architecture associated with a sexually selected syndrome in common wall lizards, by capitalizing on the adaptive introgression of coloration and morphology into a distantly related lineage. Consistent with the hypothesis that the evolution of phenotypic syndromes in vertebrates is facilitated by developmental linkage through neural crest cells, most of the genes associated with the syndrome are involved in neural crest cell regulation. A major locus was a ~400-kb region, characterized by standing structural genetic variation and previously implied in the evolutionary innovation of coloration and beak size in birds. We conclude that features of the developmental and genetic architecture contribute to maintaining trait integration, facilitating the extensive and rapid introgressive spread of suites of sexually selected characters.

GbyE: an integrated tool for genome widely association study and genome selection based on genetic by environmental interaction.

BACKGROUND: The growth and development of organism were dependent on the effect of genetic, environment, and their interaction. In recent decades, lots of candidate additive genetic markers and genes had been detected by using genome-widely association study (GWAS). However, restricted to computing power and practical tool, the interactive effect of markers and genes were not revealed clearly. And utilization of these interactive markers is difficult in the breeding and prediction, such as genome selection (GS). RESULTS: Through the Power-FDR curve, the GbyE algorithm can detect more significant genetic loci at different levels of genetic correlation and heritability, especially at low heritability levels. The additive effect of GbyE exhibits high significance on certain chromosomes, while the interactive effect detects more significant sites on other chromosomes, which were not detected in the first two parts. In prediction accuracy testing, in most cases of heritability and genetic correlation, the majority of prediction accuracy of GbyE is significantly higher than that of the mean method, regardless of whether the rrBLUP model or BGLR model is used for statistics. The GbyE algorithm improves the prediction accuracy of the three Bayesian models BRR, BayesA, and BayesLASSO using information from genetic by environmental interaction (G × E) and increases the prediction accuracy by 9.4%, 9.1%, and 11%, respectively, relative to the Mean value method. The GbyE algorithm is significantly superior to the mean method in the absence of a single environment, regardless of the combination of heritability and genetic correlation, especially in the case of high genetic correlation and heritability. CONCLUSIONS: Therefore, this study constructed a new genotype design model program (GbyE) for GWAS and GS using Kronecker product. which was able to clearly estimate the additive and interactive effects separately. The results showed that GbyE can provide higher statistical power for the GWAS and more prediction accuracy of the GS models. In addition, GbyE gives varying degrees of improvement of prediction accuracy in three Bayesian models (BRR, BayesA, and BayesCpi). Whatever the phenotype were missed in the single environment or multiple environments, the GbyE also makes better prediction for inference population set. This study helps us understand the interactive relationship between genomic and environment in the complex traits. The GbyE source code is available at the GitHub website ( https://github.com/liu-xinrui/GbyE ).

Positive selection underlies repeated knockout of ORF8 in SARS-CoV-2 evolution.

Knockout of the ORF8 protein has repeatedly spread through the global viral population during SARS-CoV-2 evolution. Here we use both regional and global pathogen sequencing to explore the selection pressures underlying its loss. In Washington State, we identified transmission clusters with ORF8 knockout throughout SARS-CoV-2 evolution, not just on novel, high fitness viral backbones. Indeed, ORF8 is truncated more frequently and knockouts circulate for longer than for any other gene. Using a global phylogeny, we find evidence of positive selection to explain this phenomenon: nonsense mutations resulting in shortened protein products occur more frequently and are associated with faster clade growth rates than synonymous mutations in ORF8. Loss of ORF8 is also associated with reduced clinical severity, highlighting the diverse clinical impacts of SARS-CoV-2 evolution.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing.

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

A generic approach to infer community-level fitness of microbial genes.

The gene content in a metagenomic pool defines the function potential of a microbial community. Natural selection, operating on the level of genomes or genes, shapes the evolution of community functions by enriching some genes while depriving the others. Despite the importance of microbiomes in the environment and health, a general metric to evaluate the community-wide fitness of microbial genes remains lacking. In this work, we adapt the classic neutral model of species and use it to predict how the abundances of different genes will be shaped by selection, regardless of at which level the selection acts. We establish a simple metric that quantitatively infers the average survival capability of each gene in a microbiome. We then experimentally validate the predictions using synthetic communities of barcoded Escherichia coli strains undergoing neutral assembly and competition. We further show that this approach can be applied to publicly available metagenomic datasets to gain insights into the environment-function interplay of natural microbiomes.

Background Selection From Unlinked Sites Causes Nonindependent Evolution of Deleterious Mutations.

Background selection describes the reduction in neutral diversity caused by selection against deleterious alleles at other loci. It is typically assumed that the purging of deleterious alleles affects linked neutral variants, and indeed simulations typically only treat a genomic window. However, background selection at unlinked loci also depresses neutral diversity. In agreement with previous analytical approximations, in our simulations of a human-like genome with a realistically high genome-wide deleterious mutation rate, the effects of unlinked background selection exceed those of linked background selection. Background selection reduces neutral genetic diversity by a factor that is independent of census population size. Outside of genic regions, the strength of background selection increases with the mean selection coefficient, contradicting the linked theory but in agreement with the unlinked theory. Neutral diversity within genic regions is fairly independent of the strength of selection. Deleterious genetic load among haploid individuals is underdispersed, indicating nonindependent evolution of deleterious mutations. Empirical evidence for underdispersion was previously interpreted as evidence for global epistasis, but we recover it from a non-epistatic model.

A hypothesis of teleological evolution, via endogenous acetylcholine, nitric oxide, and calmodulin pathways.

The Extended Evolutionary Synthesis (EES) addresses the issues in evolutionary biology which cannot be explained by neo-Darwinian theory. The EES paradigm recognises teleology and agency in living systems, and identifies that organisms can directly affect their evolutionary trajectory in a goal-directed manner, yet the physiological pathways via which this occurs remain unidentified. Here, I propose a physiological pathway via which organisms can alter their genotype and phenotype by making behavioural decisions with respect their activity levels, partitioning of resources either toward growth, defence against disease, or their behavioural response to stressors. Specifically, I hypothesize that agential, teleological decisions mediated by acetylcholine result in induced nitric oxide (NO) activity, which regulates metabolism, blood flow, and immune response. Nitric oxide, however, is also a key epigenetic molecule, being involved in DNA acetylation, methylation, and de-methylation. Further, NO alters the histone complexes which scaffold nuclear DNA strands, and is thus a good candidate in identifying a system which allows an organisms to make teleological genetic changes. The proposed mechanisms of inheritance of these genetic changes is via the paternal line, whereby epigenetic changes in the somatic Sertoli cells in animals are transcribed by mRNA and included in the germline cells - the male gametes. The microsporangium in plants, and the sporophore cells in fungi, meanwhile, are proposed to form similar systems in response to sensory detection of stressors. Whilst the hypothesis is presented as a simplified model for future testing, it opens new avenues for study in evolutionary biology.

Genomic selection in plant breeding: Key factors shaping two decades of progress.

Genomic selection, the application of genomic prediction (GP) models to select candidate individuals, has significantly advanced in the past two decades, effectively accelerating genetic gains in plant breeding. This article provides a holistic overview of key factors that have influenced GP in plant breeding during this period. We delved into the pivotal roles of training population size and genetic diversity, and their relationship with the breeding population, in determining GP accuracy. Special emphasis was placed on optimizing training population size. We explored its benefits and the associated diminishing returns beyond an optimum size. This was done while considering the balance between resource allocation and maximizing prediction accuracy through current optimization algorithms. The density and distribution of single-nucleotide polymorphisms, level of linkage disequilibrium, genetic complexity, trait heritability, statistical machine-learning methods, and non-additive effects are the other vital factors. Using wheat, maize, and potato as examples, we summarize the effect of these factors on the accuracy of GP for various traits. The search for high accuracy in GP-theoretically reaching one when using the Pearson's correlation as a metric-is an active research area as yet far from optimal for various traits. We hypothesize that with ultra-high sizes of genotypic and phenotypic datasets, effective training population optimization methods and support from other omics approaches (transcriptomics, metabolomics and proteomics) coupled with deep-learning algorithms could overcome the boundaries of current limitations to achieve the highest possible prediction accuracy, making genomic selection an effective tool in plant breeding.

A quantitative genetic model of background selection in humans.

Across the human genome, there are large-scale fluctuations in genetic diversity caused by the indirect effects of selection. This "linked selection signal" reflects the impact of selection according to the physical placement of functional regions and recombination rates along chromosomes. Previous work has shown that purifying selection acting against the steady influx of new deleterious mutations at functional portions of the genome shapes patterns of genomic variation. To date, statistical efforts to estimate purifying selection parameters from linked selection models have relied on classic Background Selection theory, which is only applicable when new mutations are so deleterious that they cannot fix in the population. Here, we develop a statistical method based on a quantitative genetics view of linked selection, that models how polygenic additive fitness variance distributed along the genome increases the rate of stochastic allele frequency change. By jointly predicting the equilibrium fitness variance and substitution rate due to both strong and weakly deleterious mutations, we estimate the distribution of fitness effects (DFE) and mutation rate across three geographically distinct human samples. While our model can accommodate weaker selection, we find evidence of strong selection operating similarly across all human samples. Although our quantitative genetic model of linked selection fits better than previous models, substitution rates of the most constrained sites disagree with observed divergence levels. We find that a model incorporating selective interference better predicts observed divergence in conserved regions, but overall our results suggest uncertainty remains about the processes generating fitness variation in humans.

Amplifiers of selection for the Moran process with both Birth-death and death-Birth updating.

Populations evolve by accumulating advantageous mutations. Every population has some spatial structure that can be modeled by an underlying network. The network then influences the probability that new advantageous mutations fixate. Amplifiers of selection are networks that increase the fixation probability of advantageous mutants, as compared to the unstructured fully-connected network. Whether or not a network is an amplifier depends on the choice of the random process that governs the evolutionary dynamics. Two popular choices are Moran process with Birth-death updating and Moran process with death-Birth updating. Interestingly, while some networks are amplifiers under Birth-death updating and other networks are amplifiers under death-Birth updating, so far no spatial structures have been found that function as an amplifier under both types of updating simultaneously. In this work, we identify networks that act as amplifiers of selection under both versions of the Moran process. The amplifiers are robust, modular, and increase fixation probability for any mutant fitness advantage in a range r ∈ (1, 1.2). To complement this positive result, we also prove that for certain quantities closely related to fixation probability, it is impossible to improve them simultaneously for both versions of the Moran process. Together, our results highlight how the two versions of the Moran process differ and what they have in common.

The cliff edge model of the evolution of schizophrenia: Mathematical, epidemiological, and genetic evidence.

How has schizophrenia, a condition that significantly reduces an individual's evolutionary fitness, remained common across generations and cultures? Numerous theories about the evolution of schizophrenia have been proposed, most of which are not consistent with modern epidemiological and genetic evidence. Here, we briefly review this evidence and explore the cliff edge model of schizophrenia. It suggests that schizophrenia is the extreme manifestation of a polygenic trait or a combination of traits that, within a normal range of variation, confer cognitive, linguistic, and/or social advantages. Only beyond a certain threshold, these traits precipitate the onset of schizophrenia and reduce fitness. We provide the first mathematical model of this qualitative concept and show that it requires only very weak positive selection of the underlying trait(s) to explain today's schizophrenia prevalence. This prediction, along with expectations about the effect size of schizophrenia risk alleles, are surprisingly well matched by empirical evidence. The cliff edge model predicts a dynamic change of selection of risk alleles, which explains the contradictory findings of evolutionary genetic studies.

The pseudoentropy of allele frequency trajectories, the persistence of variation, and the effective population size.

To concisely describe how genetic variation, at individual loci or across whole genomes, changes over time, and to follow transitory allelic changes, we introduce a quantity related to entropy, that we term pseudoentropy. This quantity emerges in a diffusion analysis of the mean time a mutation segregates in a population. For a neutral locus with an arbitrary number of alleles, the mean time of segregation is generally proportional to the pseudoentropy of initial allele frequencies. After the initial time point, pseudoentropy generally decreases, but other behaviours are possible, depending on the genetic diversity and selective forces present. For a biallelic locus, pseudoentropy and entropy coincide, but they are distinct quantities with more than two alleles. Thus for populations with multiple biallelic loci, the language of entropy suffices. Then entropy, combined across loci, serves as a concise description of genetic variation. We used individual based simulations to explore how this entropy behaves under different evolutionary scenarios. In agreement with predictions, the entropy associated with unlinked neutral loci decreases over time. However, deviations from free recombination and neutrality have clear and informative effects on the entropy's behaviour over time. Analysis of publicly available data of a natural D. melanogaster population, that had been sampled over seven years, using a sliding-window approach, yielded considerable variation in entropy trajectories of different genomic regions. These mostly follow a pattern that suggests a substantial effective population size and a limited effect of positive selection on genome-wide diversity over short time scales.

Exploring the patterns of evolution: Core thoughts and focus on the saltational model.

The Modern Synthesis, a pillar in biological thought, united Darwin's species origin concepts with Mendel's laws of character heredity, providing a comprehensive understanding of evolution within species. Highlighting phenotypic variation and natural selection, it elucidated the environment's role as a selective force, shaping populations over time. This framework integrated additional mechanisms, including genetic drift, random mutations, and gene flow, predicting their cumulative effects on microevolution and the emergence of new species. Beyond the Modern Synthesis, the Extended Evolutionary Synthesis expands perspectives by recognizing the role of developmental plasticity, non-genetic inheritance, and epigenetics. We suggest that these aspects coexist in the plant evolutionary process; in this context, we focus on the saltational model, emphasizing how saltation events, such as dichotomous saltation, chromosomal mutations, epigenetic phenomena, and polyploidy, contribute to rapid evolutionary changes. The saltational model proposes that certain evolutionary changes, such as the rise of new species, may result suddenly from single macromutations rather than from gradual changes in DNA sequences and allele frequencies within a species over time. These events, observed in domesticated and wild higher plants, provide well-defined mechanistic bases, revealing their profound impact on plant diversity and rapid evolutionary events. Notably, next-generation sequencing exposes the likely crucial role of allopolyploidy and autopolyploidy (saltational events) in generating new plant species, each characterized by distinct chromosomal complements. In conclusion, through this review, we offer a thorough exploration of the ongoing dissertation on the saltational model, elucidating its implications for our understanding of plant evolutionary processes and paving the way for continued research in this intriguing field.

Indirect genetic effects should make group size more evolvable than expected.

Group size is an important trait for many ecological and evolutionary processes. However, it is not a trait possessed by individuals but by social groups, and as many genomes contribute to group size understanding its genetic underpinnings and so predicting its evolution is a conceptual challenge. Here I suggest how group size can be modelled as a joint phenotype of multiple individuals, and so how models for evolution accounting for indirect genetic effects are essential for understanding the genetic variance of group size. This approach makes it clear that (a) group size should have a larger genetic variance than initially expected as indirect genetic effects always contribute exactly as much as direct genetic effects and (b) the response to selection of group size should be faster than expected based on direct genetic variance alone as the correlation between direct and indirect effects is always at the maximum positive limit of 1. Group size should therefore show relatively rapid evolved increases and decreases, the consequences of which and evidence for I discuss.

Dominance reversals: the resolution of genetic conflict and maintenance of genetic variation.

Beneficial reversals of dominance reduce the costs of genetic trade-offs and can enable selection to maintain genetic variation for fitness. Beneficial dominance reversals are characterized by the beneficial allele for a given context (e.g. habitat, developmental stage, trait or sex) being dominant in that context but recessive where deleterious. This context dependence at least partially mitigates the fitness consequence of heterozygotes carrying one non-beneficial allele for their context and can result in balancing selection that maintains alternative alleles. Dominance reversals are theoretically plausible and are supported by mounting empirical evidence. Here, we highlight the importance of beneficial dominance reversals as a mechanism for the mitigation of genetic conflict and review the theory and empirical evidence for them. We identify some areas in need of further research and development and outline three methods that could facilitate the identification of antagonistic genetic variation (dominance ordination, allele-specific expression and allele-specific ATAC-Seq (assay for transposase-accessible chromatin with sequencing)). There is ample scope for the development of new empirical methods as well as reanalysis of existing data through the lens of dominance reversals. A greater focus on this topic will expand our understanding of the mechanisms that resolve genetic conflict and whether they maintain genetic variation.

Landscape genomics reveals regions associated with adaptive phenotypic and genetic variation in Ethiopian indigenous chickens.

Climate change is a threat to sustainable livestock production and livelihoods in the tropics. It has adverse impacts on feed and water availability, disease prevalence, production, environmental temperature, and biodiversity. Unravelling the drivers of local adaptation and understanding the underlying genetic variation in random mating indigenous livestock populations informs the design of genetic improvement programmes that aim to increase productivity and resilience. In the present study, we combined environmental, genomic, and phenotypic information of Ethiopian indigenous chickens to investigate their environmental adaptability. Through a hybrid sampling strategy, we captured wide biological and ecological variabilities across the country. Our environmental dataset comprised mean values of 34 climatic, vegetation and soil variables collected over a thirty-year period for 260 geolocations. Our biological dataset included whole genome sequences and quantitative measurements (on eight traits) from 513 individuals, representing 26 chicken populations spread along 4 elevational gradients (6-7 populations per gradient). We performed signatures of selection analyses ([Formula: see text] and XP-EHH) to detect footprints of natural selection, and redundancy analyses (RDA) to determine genotype-environment and genotype-phenotype-associations. RDA identified 1909 outlier SNPs linked with six environmental predictors, which have the highest contributions as ecological drivers of adaptive phenotypic variation. The same method detected 2430 outlier SNPs that are associated with five traits. A large overlap has been observed between signatures of selection identified by[Formula: see text]and XP-EHH showing that both methods target similar selective sweep regions. Average genetic differences measured by [Formula: see text] are low between gradients, but XP-EHH signals are the strongest between agroecologies. Genes in the calcium signalling pathway, those associated with the hypoxia-inducible factor (HIF) transcription factors, and sports performance (GALNTL6) are under selection in high-altitude populations. Our study underscores the relevance of landscape genomics as a powerful interdisciplinary approach to dissect adaptive phenotypic and genetic variation in random mating indigenous livestock populations.